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МНОГОУРОВНЕВАЯ МОДЕЛЬ АНАЛИЗА ПРОЧНОСТИ КОНСТРУКЦИЙ 
ИЗ ПОЛИМЕРНЫХ КОМПОЗИТОВ ПРИ МНОГОФАКТОРНОМ НАГРУЖЕНИИ 


О.Г. ОСЯЕВ, Ю.А. ТАТУРИН, А.М. КОСТИН, А.В. ЖУКОВ 
(Ростовский военный институт ракетных войск) 


Предложена математическая модель для анализа сложного напряженно-деформированног О СОСТОЯНИЯ МНО- 
ГОСЛОЙНыХ несущих конструкций летательных аппаратов, представленная в виде системы уравнений СВЯЗИ 
параметров нагружения, напряжений И деформаций, Для решения системы уравнений использован кКом- 
плексный мног суровневый метод решения, ПОЗВОЛЯЮЩИЙ учесть особенности г еометрии конструкции и по- 
ведения конструкционных материалов в процессе мног офак торного нагружения. 

Ключевые слова: длительная прочность, напряженно-деформированное СОСТОЯНИЕ, ПОЛЗУЧЕСТЬ. 


Введение. На этапах эксплуатации конструкции летательных аппаратов испытывают комплекс- 
ное воздействие факторов термосилового нагружения. Особенность механического поведения 
конструкционных полимерных композитных материалов состоит в том, что при нормальной тем- 
пературе эксплуатации и сравнительно небольших уровнях напряжений они проявляют свойство 
ползучести. Поэтому при выполнении прочностных расчетов необходимо учитывать их указанные 
особенности. 

Установлено [1], что закономерности ползучести основных конструкционных полимерных 

материалов в широком диапазоне напряжений удовлетворительно описываются линейными на- 
следственными уравнениями. В общем случае пространственного теплового и напряженно- 
деформированного состояния краевая задача линейной наследственной теории ползучести сво- 
дится к решению уравнений наследственной термовязкоупругости. 
Определение напряженно-деформированного состояния конструкций из полимерных 
композитов при многофакторном нагружении. На основании принципа соответствия реше- 
ние задачи наследственной ползучести может быть приведено к решению соответствующей зада- 
чи упругости путем замены основных констант вязкоупругости материала соответствующими вре- 
менными операторами с помощью прямого символического метода Вольтерра, либо применением 
преобразования Лапласа или Лапласа — Карсона к наследственным интегралам ползучести. Для 
этой же цели могут быть использованы методы преобразования Фурье или разложения в ряды 
других видов. Выбор метода зависит от физической модели рассматриваемых процессов и свойств 
материала конструкций. 

Обозначим параметры пространственного теплового и напряженно-деформированного 
состояния (НДС), к которым применены указанные виды преобразований, верхним индексом «*», 
а временные операторы вязкоупругости, замещающие соответствующие константы вязкоупруго- 
сти материалов, — верхней чертой. Представим воздействие внешних статических, динамических 


сил и тепловых источников в виде вектора полей температур Т*, тепловых потоков 4“, напря- 


жений с” и перемещений и” на внутренней и наружной поверхностях оболочки корпуса несущей 
конструкции летательного аппарата: 


зи 


+ = ры Еы + = = = 
с ={Т Ч бен |. (1) 


Рассмотрим многослойную оболочку корпуса несущей конструкции летательного аппарата 
из вязкоупругих анизотропных материалов, отнесенную к криволинейной ортогональной системе 


координат х,х,,х,. Здесь, наряду с общей системой координат х,,х,,х,, введена также и ло- 
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кальная система х,,,х,,,х,,, связанная с каждым слоем А. При этом каждая из поверхностей 


слоя К х,,=0 совпадает с его срединной поверхностью радиусом А;’. Для каждого слоя К 
многослойной оболочки справедлива система уравнений, которая получается на основе классиче- 
ской системы нелинейных уравнений теории упругости [2] в предположении геометрической ли- 
нейности деформаций =, относительно перемещений и, и физической нелинейности материала, 


обусловленной реологическими свойствами полимерных композитов. С учетом рассмотренных 
преобразований, исходную систему уравнений представим в тензорном виде: 


до, = ‚ ди, 
о. ве И о, (2) 
дх, - . о 0%, 0% 


1 





где Х, - вектор объемных сил; Е - операторный модуль вязкоупругости; е, — деформация. 


Для перехода от уравнений равновесия к уравнениям движения компоненты вектора объ- 
емных сил представим выражением: 





2 
р Е (3) 
0, 


Применяя преобразование Лапласа к уравнению (3) и принимая параметр преобразова- 
ния р = Л ‚гдех, - эмпирическая константа материала, 


Х,=Х,-Рр 


при начальных условиях 
ди _ ди 


== (4) 


> 
И 
= 


1=0 








1=0 
получаем: 


ХХ! =Х! р! ^^" Ам Е 


‚1=1 2,3. 5 
1=0 О | (5) 





Учитывая особенности геометрии формы конструкций оболочек с использованием коэф- 
фициентов Ламе Н,,Н,,Н,, исходную систему уравнений, справедливую для каждого слоя К 


многослойной оболочки, с учетом выражений (3)-(5), запишем в следующем виде: 
— уравнения движения 








д +. ОН; + д . ОН, _. Г + 
—(Н,6)- 26, + —(Нб,)+-— а, + —(Н.Н,6\.)+ 
ах, ( 2 и) ах, 22 дх, ( 1 12) дх, 12 дх, ( тие 13) 
ее. Аи, — Ам, к =0, 
дх, #=0 0 |1=0 
д . ОН ры 20 . ОН, _, Г А 
— (Но. - 6 +—(Н,в..) +в, + —(Н,Н,6,.)+ 
дх, ( 1 >) дх, 11 ах, ( р) 12) дх, 12 дх, ( 12, эз) 
6) 
ОН, + . . ди ( 
+Н, 03+) -р| ^^, Ам, - —— =0, 
. =0 0 |1=0 
5 (Майи) = © + == О) $ -—(Н.Н,озз) са 
ья дх. дх› дх. 
НГ * * * д 
ааа р Л? — Ам, а =0; 
ь (=0 ОЕ |1=0 
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— физические уравнения связи напряжений и деформаций 











ож Ур Уз * , 1, 
5 — © —@ —^о &, =— 0 
11 И 22 337 12 12 7 
Ей Е» Езз С 
р о = Мы 4. У - * 1, 
= 12 23 =: 
8> =: Е 9 Е О Е Оз 7 83 — с Оз 7 (7) 
22 22 33 13 
- > „ У * 1 
13 23 ы * 
5... =—— О — © =—-20. $. = —_б 
33 33 И 221 23 23 
Езз Езз Езз С»; 


Здесь характеристики материала Е и С представляют операторные модули вязкоупруго- 
сти и вязкого сдвига соответственно, у-— коэффициент Пуассона; 
— геометрические уравнения Коши 
‚1 бы 1 он. 1 ОН. 
&и = + 2+ и, 
Н.Н, дх, 





= и 
Во ЗН 
о, на ЗО 
5>> = Е 5+ и, 

Н,Н, 9х, 





= и 
Н, @, Н,Н, дх, 
мы Гон. тан, . 
- + 























8: = и и, (8) 
Н, Хх. НН, 9% НН. 9х, 
Ро Е ВИ Е 
52 = + ; 
Н, 9х, Н, 9х, 
е.. Тоби: 0: 
&3 = + , 
Н, @х, Н, 9х 
. д. Тб 
83 = + , 
Н, 9х, Н, 0х, 
— начальные условия 
жа 65`\° |(05 ре 
сб =5 а ба | =рб -б ; (9) 
1=0 0 ОЕ =0 =0 
— граничные условия 
на наружной и внутренней поверхностях пакета слоев: 
* + * * + * * + * 
О = (Оз) , 91 = (65:1) , 93 = (зв) , 
(10) 
бк = (бв»), бк =(б,), Обр =(би); 
поверхностях контакта смежных слоев: 
И Иды, Оззд = Оззи 
(11) 
бук Оу, 51 =1 2,3, Е=1, 2, ..,К; 
боковых поверхностях смежных слоев: 
* * * 1 * * А 
О (+ енк) +9124 [5 № Эз 2Оид, 
* 1 * * * * Го 
О оз Е Э:з* Орд (+ >> + )= Ол, (12) 


* * 


Од оз Эрл | ТОрь ое Уи | ТОбзх =Овд, 


9724 ое = Чик | Обрх оз +9 ТО Ок. 
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Для решения системы уравнений (6)-(12) выберем в качестве неизвестных функции, с по- 
мощью которых выражаются условия контакта между слоями. 

Компоненты тензоров деформаций, а также слагаемые, содержащие величины углов по- 
ворота в уравнениях (12), определяются с помощью известных соотношений нелинейной теории 
упругости [2]. Решая систему (6)-(8) относительно этих функций, получаем для ортотропной ци- 
линдрической оболочки следующую разрешающую систему уравнений: 

. . а 2 
- = то; А 7 : ы т +9. )- те ыы (Аз +0, )- 

















бз 11 вех 
1х, х Хх х 1х. 
0х, 3 ° @% з 0% 
[9 . . ОЕ ‚ ОГ . 
трх. Е (Аи -А; а» )- — Аз ие: -Х, 
х. дх 0х, | | о О = ” дх 
. ди 
+Р| Ам, — Ам, о , 
=0 ог. [0 
до’. а ое + 
ВО = 45 а (А +А.,)- 
дх, %: х. дх. х, дх 0х, 
ОА Ах 1 . 
2 2 
= Я? (Аз +9. )- 2 (А; +9, )+ 20 22.012 = 
1 хз ОХ, 3 
1 ди". Вы ЭР $ 
= а (АЗ т 20 >> о )- — 4, „› их Х, 
Хз ОХ, Хз Хх? 
. ди 
+Р| м. — Ам, -—^ ; (13) 
=0 ОЕ 1=0 
до, до, 1 [ дс, . $ 95 ‚ би . 
== — т Сзз А, 22033 -А,,>> == А, Г |+ 
дх, Об | | : | 
Же а и. = 0?и. . 
2 3 3 
Ух (^; „> +09>20 +4); „зи; ) 10 2 22,0 2 Х, 
х; дх, 9х! дх> 
. ди 
+р| Ам; — Ам, - 
в=0 Ге] 





* 


ОИ ди: ЗОВ. О, жж ь (1 


а ЕЛЕ 
' 44 23 
[бя [62% дх, 


и. * * * и” * * и * 

ди, ди 1 ди, 

= 71 33063 + И, 33 + — И зз| + +т, Г 

о ее д 
Напряжения с',, с',, о’, определяются с помощью физических и геометрических урав- 

нений (7), (8): 


* 


. . вый . 
ОА | р 
Хх Хх 


* 


— 933 (аи ча ) +4 иГ, 


2,11 3 
2 


| 3 
ба = А [+ ты + 0 р (14) 
Же и 


. и, аа ‚ ди, (++ а . 
= 1 2 
о». =Д,>> р. +— А, | и + Ре В бзз (а.А! > +аА )+ А, »Т. 


2,22 2,22 
т № 2 
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Коэффициенты в уравнениях (13) и (14) определяются следующими соотношениями меж- 
ду физико-механическими и теплофизическими параметрами слоев оболочки корпуса летательно- 
го аппарата: 


* * * * * * * * * * 2 * 
о >> Ч. же Я! Ч, ее Я! Ч, а а @ бб. .ы ана» — (а>)_. 
мя г Аи Яя г А12 = ор А22 = г Аза = х 
А А А А А 
* * * * * * * * *ж * * * * * * * 
А =ав [а — (а,›) ] ‚ Аи= -а А, и-@3А и; А» =-азА»›-азА.»; Аи = -9А, и; 
* * * * * * * * * * * * 
А» =0 А —0.Д, > ’ Пзз =азА а: +а3Д, > , Й, 3 =азА, 1 +азА. 2; (15) 
* * * * * * : * * * Е 
т = ЧП —аз А, +аз, Т› 3 = И, 33 — 952533 + Оз, 
+ 1 й 1 + 1 * У ы У * У 
—_ : № } № 2 Аи Е УЕ 
а ==, “2 = ‚ Чжу @= т’ аз = 2. аз = Е. ^ 
и 22 33 22 33 33 
* 1 ы * 1 к * 1 
о 4-6. №6. 
23 13 12 


Представим неизвестные функции б”, а также параметры нагрузки б* в виде двойных 
тригонометрических рядов: 


Се] [е] 
[м аи = а те Г х с0$их,, 
=1 и=0 


С.) со 
* * + ОВ * * Е . тт . 16 
{и 0»0} = о неба т Г Ж Ш Их, , (16) 
о 


т=Ё п= 


З 


[3 + 
т [М 


* * + И * * ен . тд 5 
[и › 913505; г. {и тп > О 33,ти ? О 3Зти т 7 ж с05 Их, , 


Я 
Ь 


п= 


2п 


1 
СЕ 
913” = о сот С05 пх. фах, Г. 
оо 
12п 
+ сы Е 
б>3„„ = [[с 03 их ших, фах. , (17) 
00 
12 
4 
933. ии =| 63; $ тих с05 их ах ах, . 
00 


Производные по времени разложим в конечные разности с шагом Л: 


* 





05” _25'(1)-50 (1.)+465'(&_,)-5’(4,.) 
д ДЕР Г 
(18) 





05’ _ 36 (1)-45'(1)+5 (1...) 

и 2 
Подставим разложенные в ряды (16) и конечные разности (18) искомые функции 5’ в 
разрешающую систему уравнений (13). Получим систему обыкновенных дифференциальных 
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уравнений для каждого временного шага 1, и каждой пары волновых чисел т и п разложения в 


двойные ряды Фурье соответственно по продольной х и окружной х, координатам: 
до" (1.) 6, (1.) №: д * * ти 
а =. — АА „бзз (1, + и, (Аи + 3,12 р био + 
Ох ра а 
п? > > 
ео) — Ям, (1, У А, И +А, а — Я ыиз (1) — 4 1„Га,)- 


3 Хх Хх 


-.(,)+- ти, (11) -5и Г +4 (1(_,)- и! (3)|; 








до: (1) 25) в п + + й 
= + А,2>633(1,) - ГА „и, (1, э(^, 1,22 + А; 12 ) + 
[бя х; х; х. 

, еси ей 1 
чи, ((,)(^; и А, +01, А» + —0>2.0 +5020 ) + 
а. № Ж 


п * ИП + 
Е. +260 ) г 8) — (19) 
3 


3 


—х, (5. ее дети, = 5и> (1, _ ‚+4. (1,_ = и (1 ЗИ 








00-3 (1 0,3 (1 03), Де 34а 
до. (1,) )_ ое ) ы И ие (4 19) 1,22 Ани, (#)+ 
дх. ж х; а 
п. + й 1). 1 .| Та.) 
+ и, (.) (4) „> ы О + Е А. бои +5 боди" |+ 
№ № № № 


-Х; (1, ея диз, )- 5из (Е +45 (Е 2) ве. $]; 












Е 2 ра Я 44 ©) (1, в и> (1, уе и (1, ); 
ух т № 
я () * * №: * И * * 1 * * * 
= 33033 (,)- ПМ (+ П, 33 = Й) 33 (,)+ т, 3[ (,). 
&%, № № 


Решение системы уравнений (19) осуществляется численным методом дискретной ортого- 
нализации [3], который позволяет автоматически удовлетворять граничным условиям контакта 
слоев оболочки. Статическое напряженно-деформированное состояние оболочки, обусловленное 
постоянными или медленно изменяющимися нагрузками в процессе эксплуатации, определяется с 
помощью методов расчета наследственной ползучести [1]. При воздействии динамических нагру- 
зок на рассматриваемую оболочку корпуса, полученные значения с, и с„, ‚, характеризующие 
предыдущее статическое состояние конструкции, используются в качестве параметров предвари- 
тельного нагружения в уравнениях (19). Результатом решения системы уравнений (19) являются 


функции с,„. Искомые функции с’, определяют с помощью выражений (16) путем двойного 
суммирования с’. „. Оставшиеся искомые напряжения с\,, о',, о», находят также путем двойно- 


го суммирования результатов разложения этих функций, определяемых уравнениями (14), анало- 
гично уравнениям (16): 


Аз ны в 
си = Уд 1, О и (и; тяП + Из >" А, АИ тя 5тА„х, со$их, , 
т=1 п=0 № 


© 05 1 
* № ай * к . 
012 = уз И ттт -. —_ ИП с037.„х, Ш пх, ! (20) 
т=Ё п=0 Х, 
* 
* * * 2,22 * * * ча * . А” 
= у о 2" (и „п Йа ) — А, АИ т | А, их, С0$ ИХ, . 
т=Ё п=0 3 
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Блок-схема алгоритма решения задачи приведена на рисунке. На первом шаге осуществ- 
ляется ввод данных о массово-геометрических, физико-механических, теплофизических и дефор- 
мационных характеристиках материалов и конструкции корпуса летательного аппарата. На вто- 
ром шаге определяются коэффициенты, учитывающие данные свойства материалов и конструк- 
ции, с помощью соотношений (15). Одновременно на шаге 3 выполняется операция задания на- 
грузки, действующей на конструкцию, выраженной через наружные и внутренние напряжения, 
перемещения и температуру. Заданные составляющие функций параметров нагружения исполь- 
зуются для определения коэффициентов пошагового разложения нагрузки по волновым числам т 
и л в меридиональном и окружном направлениях соответственно. Для этого используются двой- 
ные интегральные выражения для коэффициентов анализа Фурье (17). Далее, на шаге 5 с помо- 
щью выражений (16) определяются функции разложения нагрузки и параметров напряженно- 
деформированного состояния в двойные тригонометрические ряды. 

















1. Задание 2. Определение 
параметров коэффициентов 


конструкции 
Ель Е›> Езз О А;, Ау, 
Св 23 ту, Пу 


У12, У13, У23 @1Ь @2> 
@33 рь 
$ Е _› В, И 


3. Задание 4. Определение 5. Определение параметров 


нагрузки функций НДС для пошагового 


т т =: ы т т интегрирования 
= у в г ва , 13, ти’ © 23ти › 9 33. ти › 
м + + + б1з,лт’ 9 23.ит › 9 ЗЗ,ти* 
2ы в 33 и В и й и а 
И, Но, И, Ти ти? "2. тп” " З.ти р 

























6. Задание 7. Определение параметров 8. Определение 
НДС 
о д бб ба, 9. и 
911, ти’ 9 22 тт у 9 33,ти’ 
То, то, по 923, б 12, И. И2 Из 
Нь, Нь б12,ти’ 9 1Зтт” 9 23 ,ти ' 


Иттю И2ть ИЗть 














10. Вывод 
результатов 
на дисплей 
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Физико-математические науки 








На шестом шаге алгоритма задаются параметры численного интегрирования (начальные 
значения температуры, волновые числа, число шагов интегрирования по толщине каждого слоя и 
число точек дискретной ортогонализации). На шаге 7 выполняется решение системы обыкновен- 
ных дифференциальных уравнений (19). На шаге 8 определяются искомые полные функции на- 
пряжений и перемещений путем двойного суммирования волновых функций (20). Для решения 
уравнений (19) и (20) используются результаты, полученные на шагах 2, 5 и 6 алгоритма, как по- 
казано на блок-схеме (см. рисунок). 

Системы уравнений (6)-(20) и рассмотренные методы их решения могут быть использова- 
ны для решения физически нелинейных задач термоупругости и термовязкоупругости при много- 
факторном статическом и динамическом нагружении многослойных оболочек сложных несущих 
конструкций летательных аппаратов из полимерных композитов в условиях эксплуатации. 
Заключение. Представлена математическая модель для анализа сложного напряженно- 
деформированного состояния многослойных несущих конструкций летательных аппаратов. Реше- 
ние системы уравнений осуществлено комплексным многоуровневым методом, позволяющим 
учесть особенности геометрии конструкции и поведения конструкционных материалов при мно- 
гофакторном нагружении. 
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